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Abstract 

Numerical simulations are used to determine the influence of the 
non-local and local interactions on the intermittency corrections in 
the scaling properties of 3D turbulence. We show that neglect of 
local interactions leads to an enhanced small-scale energy spectrum 
and to a significantly larger number of very intense vortices ( "torna- 
does") and stronger intermittency (e.g. wider tails in the probability 
distribution functions of velocity increments and greater anomalous 
corrections). On the other hand, neglect of the non-local interactions 
results in even stronger small-scale spectrum but significantly weaker 
intermittency. Thus, the amount of intermittency is not determined 
just by the mean intensity of the small scales, but it is non-trivially 
shaped by the nature of the scale interactions. Namely, the role of 
the non-local interactions is to generate intense vortices responsible 
for intermittency and the role of the local interactions is to dissipate 



them. Based on these observations, a new model of turbulence is pro- 
posed, in which non-local (RDT-like) interactions couple large and 
small scale via a multiplicative process with additive noise and the 
local interactions are modeled by a turbulent viscosity. This model 
is used to derive a simple toy version of the Langevin equations for 
small-scale velocity increments. A Gaussian approximation for the 
large scale fields yields the Fokker-Planck equation for the probability 
distribution function of the velocity increments. Steady state solu- 
tions of this equation allows to qualitatively explain the anomalous 
corrections and the skewness generation along scale. A crucial role is 
played by the correlation between the additive and the multiplicative 
(large-scale) process, featuring the correlation between the stretching 
and the vorticity. 
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A puzzling feature of three-dimensional turbulence are the large devia- 
tions from Gaussianity observed as one probes smaller and smaller scales. 
These deviations are usually believed to be associated with the spatial in- 
termittency of small-scale structures, organized into very thin and elongated 



intense vortices ( "tornadoes" ) 0, ^ . They are responsible for anomalous 
corrections to the normal scaling behavior of structure functions associated 
with the Kolmogorov 1941 |25] picture of turbulence. In this picture, energy 
containing structures (the so-called "eddies") at a given scale, interact with 
other eddies of smaller but comparable size to transfer energy at a constant 
rate down to the dissipative scale. A simple prediction of this local picture 
of turbulence is the famous k~^/^ energy spectrum, which has been observed 
in numerous high Reynolds number experimental data and numerical simu- 
lations. The local theory of turbulence further leads to the prediction that 
the n^^ moment of a velocity increment Su^ = u{x + i) —u{x) over a distance 
i should scale like i"'^^. This behavior has never been observed in turbu- 
lent flows, and is it now widely admitted that anomalous corrections exist 
at any finite Reynolds number. Various scenarii have been so far proposed 
to explain and compute the anomalous corrections. To mention but a few: 
spatial intermittency of the energy dissipation [2^, multi-fractal scaling [|C 



large deviations of multiplicative cascades pT|, extremum principle [|TT|, zero- 
modes of differential operators , scale covariance |]I5| . These approaches 
all try to model the breakdown of the exact local scale invariance underlying 
the Kolmogorov 1941 picture. In a recent study of finite size effects, DubruUe 
[H] showed that some properties of the structure functions (non-power law 
behavior, nonlinear exponents,..) could be explained within a framework in 
which finite size cut-offs play a central role, and are felt throughout the so- 
called "inertial range". Such finding is in clear contradiction with the "local" 
K41 theory, in which eddies in the inertial range are insensitive to the UV 
and IR end of the energy spectrum, and only interact with their neighbors (in 
the scale space) via "local" interactions (involving triads of comparable size). 
This observation motivated us to consider a new scenario for turbulence, in 
which anomalous corrections and deviations from Gaussianity are the result 
of non-local interactions between energy containing structures. By non-local 
interactions, we mean interactions between well separated scales (or highly 
elongated wavenumber triads). A numerical analysis on the role of the dif- 
ferent triadic interactions in the energy cascade have been previously done 
by Brasseur et al [0 and Domaradzkiet al |T3|. We have recently demon- 
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strated via high resolution numerical simulations that in 2D turbulence, the 
small scale dynamics is essentially governed by their non-local interactions 
with the large scales |]^. This feature seems natural in view of the 
large-scale condensation of vortices 0. As a result, the weak small scales are 
more influenced by the strong large-scale advection and shearing, than by 
mutual interactions between themselves. This makes non-local interactions 
the dominant process at small scales. 

2D turbulence is very special, because there is no vortex stretching. As a 
result, there is no increase of vorticity towards smaller scales as is observed 
in 3D turbulence. A quantitative feature underlying this difference is the 
shape of the energy spectra at small scales: it is in 2D turbulence, much 
steeper than k~^^^ 3D energy spectra. In fact, a simple estimate shows that 
the borderline between local and non-local behaviors is precisely this 
spectrum: only for energy spectra steeper than can one prove rigorously 
than the dominant interactions are non-local. It is therefore not our intention 
to claim that 3D turbulence is non-local, and in fact we do believe that local 
interactions are responsible for the k~^^^ energy spectra. However, it is not 
unreasonable to think that evolution of the higher cumulants (responsible for 
deviations from Gaussianity at large deviations) is more non-local than it is 
for the energy spectrum. Indeed, calculation of higher cumulants in Fourier 
space involves integrations over larger sets of wavenumbers which (even if 
close in values pairwise) cover larger range of scales than in the case of lower 
order cumulants. 

A natural tool to study the role of the non-local interactions is the numer- 
ical simulation, because it allows a direct check of their influence by artificial 
switching-off of the elongated wavenumber triads in the Navier-Stokes equa- 
tions or, on the contrary, retaining only such triads. A limitation of this 
approach lies in the restricted range of Reynolds number we are able to sim- 
ulate. However, a recent comparison showed that anomalous corrections 
and intermittency effects are quite insensitive to the Reynolds number. This 
shows relevance of a low Reynolds number numerical study of intermittency. 
This approach is detailed in the next section, where we examine the dynam- 
ical role of the local interactions at the small scales using a simulation in 
which these interactions have been removed. We show that, as compared 
with a full simulation of the Navier-Stokes equations, such "nonlocal" simu- 
lation exhibits a flatter spectrum at small scales and stronger intermittency. 
As a qualitative indicator of intermittency we use plots of the vortex struc- 
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tures whose intensity greatly exceeds the r.m.s. vorticity value ("tornadoes") 
and PDF's plots, whereas to quantify intermittency we measure the structure 
function scaling exponents. We show that main effect of the local interactions 
can be approximately described by a turbulent viscosity, while the non-local 
interactions are responsible for development of the localized intense vortices 
and the deviations of Gaussianity. To validate that the enhanced intermit- 
tency is not merely a result of the increased mean small-scale intensity (which 
could be also caused by some reason other than the non-locality) we perform 
another numerical experiment in with the elongated triads were removed. 
Such "local" experiment resulted in even stronger small-scale intensity but 
the intermittency became clearly weaker. In the Section 2, we show how our 
findings can be used to provide both a qualitative estimate of the intermit- 
tency exponents, and a derivation of the Langevin equations for the velocity 
increments. 



1 Non-local interactions in 3D turbulence 

1.1 The problematics 

We consider the Navier-Stokes equations: 

^tu + u • Vu = -Vp + i/Au + f , (1) 

where u is the velocity, p is the pressure, u is the molecular viscosity and f 

is the forcing. In a typical situation, the forcing is provided by some bound- 
ary conditions (experiments) or externally fixed, e.g. by keeping a fixed 
low-wavenumber Fourier mode at a constant amplitude (numerical simula- 
tion). This situation typically gives rise to quasi-Gaussian large-scale velocity 
fields, while smaller-scale velocities display increasingly non-Gaussian statis- 
tics. Presence of the forcing guarantees existence of a stationary steady state 
in which the total energy is constant. In absence of forcing, the turbulence 
energy decays steadily, due to losses through viscous effects. However, start- 
ing from a quasi-Gaussian large-scale field, one can still observe development 
of increasingly non-Gaussian small scales in the early stage of the decay. We 
will study the effect of the non-local interactions on the statistics of such 
non-Gaussian small scales. For this, we introduce a filter function G'(x) in 
order to separate the large and small scales of the fiow. In our numerical 



5 



procedure, the filter G will be taken as a cut-off. We have checked that the 
results are insensitive to the choice of the filter, provided the latter decays 
fast enough at infinity. Using the filter, we decompose the velocity field into 
large scale and small scale components: 

u(x,t) = U(x,t) + u'(x,t), 

U(x,t) = M = y"G(x-x')u(x',t)rfx'. (2) 

Equations for the large scales of motion are obtained by application of the 
spatial filter (|[) to the individual terms of the basic equations (0). They are: 

d^U, = 0, 

dtUi + djUiUj + djUiUj + UjUi + djUjU] = 

-d,P + uAUi + Fi, (3) 

In these equations, we have dropped primes on sub-filter components for 
simplicity; this means that from now on, any large-scale quantities are de- 
noted by a capital letter, while the small-scale quantities are denoted by a 
lower case letter. Equation for the small-scale component is obtained by 
subtracting the large-scale equation from the basic equations (|Tp; this gives 

djUj = 0, 

dtUi + dj ((f/i + Ui){Uj + Uj) - {Ui + Ui){Uj + Uj)) = 

-dip + uAui + fi, (4) 

Several terms contribute to the interaction of scales: non-local terms, involv- 
ing the product of a large scale and a small scale component, and a local 
term, involving two small-scale components. One way to study the dynam- 
ical effect of these contributions at small scale, is to integrate numerically a 
set of two coupled equations, in which the local small scale interactions have 
been switched off at the small scales 0. This corresponds to the following set 
of equations, 

dtUi + djUiUj = —djUiUj — djUjUi — djUiU] 

^We do not switch off the local interactions at the large scales: this would hinder the 
cascade mechanism and prevent small scale generation from an initial large scale field. 
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-d,P + uAUi + U (5) 
dtUi + dj{UiUj) + dj{uiUj) = -diP + uAui + ai, (6) 
9jUj=dju, = 0, 

where 

ai = (TWj - UiUj + UjUl + U~u^ . (7) 

The later describes a forcing of the small scales by the large scales via the 
energy cascade mechanism. This term is always finite even when the external 
forcing f (which is always at large scales) is absent. The small scale equa- 
tion is linear and it resembles the equations of the Rapid Distortion Theory 
We shall therefore refer to this new model as the RDT model. The 
corresponding solution was then compared with a reference simulation per- 
formed at the same resolution, with the same initial condition. Note that this 
comparison is rather expensive numerically: splitting the equations of mo- 
tions between resolved and sub-filter component leads to additional Fourier 
transforms, and increases the computational time by a factor 3. This sets a 
practical limitation to the tests we could perform on our workstation. Also, 
an additional limitation came from the need of scale separation between the 
"large" and the "small" scales. This scale separation is mandatory in or- 
der to define "non-local" interactions. Their influence on anomalous scaling 
can be checked only if the typical small scale lies within the inertial range. 
For this reason, we were led to consider a situation of decaying (unforced) 
3D turbulence, with a fiat energy spectrum at large scale, and an "inertial 
range" mostly concentrated at small to medium scales (10 < A; < 40 for a 
256^ simulation). Indeed, at the resolution we could achieve, forced turbu- 
lence developed an inertial range of scale around = 8, too small for the 



scale separation to be effective (see p9| for a study and discussion of this 



case and its relevance to LES simulations). Decaying turbulence does not, 
by definition, achieve a statistically stationary state, with mathematically 
well defined stationary probability distribution functions (PDF). Therefore, 
all the PDF were computed at a fixed time which we have chosen to be at 
the end of each simulation (at t=0.48). 
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1.2 The numerical procedure 



1.2.1 The numerical code 

Both the Navier Stokes equation (|lj) and the set of RDT equations (||) were 
integrated with a pseudo-spectral code (see for more details on the code). 



In the RDT sharp cut-off in Fourier space was used to split the velocity 

field into large and small-scale components in Fourier space and all non- 
linear terms were computed separately in the physical space. The aliasing 
was removed by keeping only the 2/3 largest modes corresponding to the 
85 first modes in our case. The calculations presented here were done with 
256^ Fourier modes and a viscosity of 1.5 10~^ corresponding to a Reynolds 
number 57 < Rx < 80 (where Rx is the Reynolds number based on the Taylor 
micro-scale A). 



1.2.2 The simulations 

The test was performed in a situation of decaying turbulence were the forc- 
ing term /j is set to zero. The initial condition was computed by a direct 
simulation using a Gaussian velocity distribution as a starting field. We have 
chosen the velocity field after several turnover times to be the initial condi- 
tion for our simulations: the direct numerical simulation (DNS), the RDT 
simulation (eq. |^) and the "local" experiment (see the end of this section). 
In order to allow enough energy at the large scales, the sharp cut-off filter 
was taken at the wavenumber k = 24 corresponding to approximatly 5 Kol- 
mogorov scales. Because of the very low speed of the RDT simulation (three 
times more expensive than the DNS) the two simulations were performed 
between t=0 and t=0.48 corresponding to approximately 2.5 turnover times. 



1.3 Comparison of the RDT and DNS experiments 
1.3.1 Spectra 

The comparison of spectra is shown in Fig In the DNS case, one observes 
a classical evolution, in which the large scale energy decays while the inertial 
^-5/3 Yange tends to move towards smaller scales. It can be seen that the 
inertial range (characterized by the —5/3 slope) only marginally exists. In the 
RDT case, one observes a similar evolution at large scale, while a tendency 
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Figure 1: Comparison of energy spectra at t=0.48 obtained via the DNS and 
the RDT simulation. 



towards a flatter spectrum is observed near and beyond the separating scale 
(beyond which local interactions are ignored). We checked that this behavior 
is not sensitive to the resolution. 

The range of computed scales is insufficient to find a reliable value of 
the spectral slope in the RDT case. However, this slope can be predicted 
by a simple dimensional argument. This argument was presented for the 
2D case in but it is essentially the same for 3D. Indeed, in the RDT 
case the small-scale equations are linear and, therefore, the energy spectrum 
E{k) must be linearly proportional to the energy dissipation rate, e. In 
this case, the only extra dimensional parameter (in comparison with the 
local/Kolmogorov case) is the large-scale rate of strain a. There is the only 
dimensional combination of e, a and wavenumber k that has the dimension 
of E{k); this gives 

= —k-\ (8) 

a 

where C is a non-dimensional constant. In our case, our resolution is too 
low to be able to check whether the RDT spectrum follows a k~^ law, but 
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we clearly see the tendency to a flatter than —5/3 slope. The RDT case is 
reminiscent of the boundary layer, in which a spectra have been observed 
. This is not surprising because presence of the mean shear increases the 
non- locality of the scale interactions corresponding to RDT. In fact, an exact 
RDT analysis of the shear flow does predict formation of the spectrum 



m 



1.3.2 Structures 

3D turbulence is characterized by intense thin vortex filaments ( "tornadoes'" 



P, |9, 50, i3]. Their radii are of order of the dissipative (Kolmogorov) scale 
which in this case is determined by the balance of the large-scale straining 
and viscous spreading. In this respect, these vortices are similar to the clas- 
sical Burgers vortex solution. Fig. ^ shows a comparison of the "tornadoes" 
observed in the DNS and in the RDT which visualized by plotting surfaces 
of strong vorticity (|a;| > 3.5 ujrms)- In both cases, thin filamentary struc- 
tures are observed, but they appear to be much more numerous in the RDT 
case. Obviously, local interactions tend to dissipate the "tornadoes" which 
can be interpreted as a mutual distortion and entanglement of "tornadoes" 
preventing their further stretching by the large scales. On the macroscopic 
level, this can be regarded as an additional, "turbulent" , viscosity produced 
by the local interactions. This is compatible with the flattening of the en- 
ergy spectrum in the RDT case, which we interpreted above in terms of the 
turbulent viscosity effect. 

It is interesting that the Burgers vortex is essentially a linear solution 
because of the cylindrical shape of this vortex which prevents appearance of 
the quadratic (in vorticity) terms. Such a linearity is a typical feature of all 
RDT solutions. On the other hand, there is another candidate which has 
often been considered to be responsible for intermittency: this is a vortex 
reconnection process which is believed to lead to a finite time singularity 
formation (at least for inviscid fluids). Note that the vortex reconnection 
is an essentially nonlinear process in which the local scale interactions are 
playing an important role and cannot be ignored. Indeed, there is no finite 
time blow-up solutions in linear RDT. Likewise, the vorticity grows only 
exponentially in the Burgers vortex and it does not blow up in a finite time. 
From this perspective, our numerical results show that the local (vortex- 
vortex) interactions mostly lead to destruction of the intense vortices and 
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Figure 2: Comparison of vortex structures (isoplot of vorticity 
\iu\ = 3.5 ujRMs = ujc at t=0.48) for a) DNS, b) RDT, c) RDT + con- 
stant turbulent viscosity, d) RDT + the RNG-type turbulent viscosity. The 
number of "tornedoes" estimated by = pdf{\uj\) duj are respectively 
8245 (a), 21669 (b), 10226 (c) and 10779 (d). ' 

prevention of their further Burgers-like (exponential) growth which has a 
negative effect on the intermittency. This process seems to overpower the 
positive effect of the local interactions on the intermittency which is related 
to the reconnection blow-ups. At the moment, it is not possible to say if the 
same is true for much higher Reynolds number flows. 

1.3.3 Turbulent viscosity 

Comparison of the DNS and RDT results for the time evolution of the total 
energy is shown in (Fig. ^. One clearly observes a slower decline of the total 
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energy in the RDT if there was a lower viscosity. This resuhs is not 
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Figure 3: Comparison of the time evolution of total energy. 



surprising: it is well known that the influence of energy motions onto well 
separated large scales 



see 



)T| for systematic expansions) is through 
an effective eddy viscosity, supplied by the < uu> term. Our result suggests 
that, to a first approximation, the difference between the RDT and the DNS 
could be removed by including an additional "turbulent" viscosity in the 
RDT simulation. For the sake of simplicity, we decided to choose an isotropic 
tensor, chosen as to conserve the total energy. We tried two simple viscosity 
prescription: one in which Vt is constant, and one (Fig. ^ in which the 
viscosity prescription follows the shape dictated by Renormalization Group 
Theory (see e.g. ): 



q '^E{q)dq 



1/2 



(9) 



The constants were adjusted as to obtain a correct energy decay (Fig. |^). 
They are z/j = 0.0002 for the constant Ut prescription, and A = 0.02 for 
the other one. We elaborate more on the choice of this turbulent viscosity 
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Figure 4: Spectrum of the turbulent viscosity computed by at t=0.48 

in Section |2.3| . Yet another method we tried was to replace the neglected 
nonlinear term (interaction of small scales among themselves) with its mean 
value. Dividing this mean nonlinear term by fc^ for each k one can compute 
the turbulent viscosity t't(fc). The result is interesting: ut turns out to be 
nearly independent of k; this provides an extra justification for the simple 
model in which z/^ = const. 

The energy spectra and energy decay obtained with this new RDT sim- 
ulation are shown on Figs. ^ and ^ One sees that one now captures exactly 
the energy decay of the DNS. The energy spectra become closer to the DNS 
result at low and intermediate k, whereas at the high k they depart from 
DNS indicating an over-dissipation of the smallest scales. The later is an 
artifact of our crude choice for the turbulent viscosity ignoring its anisotropy 
and possibility for it to take negative values. The turbulent viscosity also 
influences the anomalous corrections, as we will now show it. 
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Figure 5: Time evolution of total energy obtained via the DNS and the RDT 
simulation with the two different turbulent viscosities. 



1.3.4 PDF's and exponents 

We conducted a statistical study of our velocity fields corresponding to the 
end of the simulations (t = 0.48). As usual for study of the anomalous 
properties of turbulence, we consider the velocity increments over a distance 
1, 

(5ui = u(x + l) -u(x). (10) 

As usual, we will deal with the longitudinal and transverse to 1 velocity 
increments, 6ui\\ = {6ui ■ 1)/1 and 6ui± = {6ui x 1)/1 respectively, where 
/ = |1|. Figs. 1^, Hand|^ show the probability distribution functions (PDF) of 



the longitudinal increments and Figs. [10], |TT| and the PDF of transverse 
increments , for 3 values of /, obtained by the DNS, and our different RDT 
simulation (with and without turbulent viscosity). 

At large scale, one observes a quasi- Gaussian behavior, with the develop- 
ment of wider tails as one goes towards smaller, inertial scales. This widening 
of the PDF's is a classical signature of the anomalous scaling observed in tur- 
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Figure 6: Comparison of energy spectra obtained via the DNS and the RDT 
simulation with two different forms of the turbulent viscosity 



bulence. It can be measured by studying the scaling properties of the velocity 
structure functions, 

5p(£) =< K > • (11) 

In the inertial range, the structure function vary like For low Reynolds 
number turbulent flows, the scaling behavior in the inertial range is very 
weak or undetectable because the inertial range is very short. To exemplify 
this point, we show in Fig. |13| and |l^ the structure functions as a function 
of the scale separation for the longitudinal and the transverse velocity incre- 
ments (the structure functions from the RDT simulation were shifted by a 
factor 10 for the clarity of the figures). Given the very weak scaling of our 
structure functions, we may use the extended self-similarity (ESS) property 
[|], which states that Sp{i) ~ 5*3'"'''^ even outside the inertial range of scales. 
We use this property because it allows to find the scaling exponents in a 
more unambiguous way |^. The measured exponents in our DNS are shown 
in Fig ^ and in Table |T[ In both the longitudinal and transverse case, they 
are in agreement with the previously reported exponents and they dis- 
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Table 1: scaling exponent of the velocity structure functions measured using 
the ESS property, in various simulations, at t=0.48. 
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RDT + Vise. 1 


RDT + Vise. 
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0.721 
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1. 


1. 
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1.843 
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1.991 
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Figure 7: Comparison of the PDF of the longitudinal velocity increments 
defined by Eq. (|10D for I = 27r/256 (5. 10^ statistics at t = 0.48 for the 
velocity field from the four simulations: DNS (circle), RDT (crosses), RDT 
+ constant viscosity (diamonds) and RDT + viscosity computed from RNG 
(triangles) 



play a clear deviation from the "non-intermittent" value (p = p/3. However 
the difference between transverse and longitudinal exponents appear to be 
somewhat larger than the one observed by Dhruva et al [0 and Camussi 
et al ip. Corresponding quantities for the RDT simulation are shown in 
Fig |T5| and Table |^. One sees that the RDT statistics display larger and 
more intermittent PDF tails for small scales, which makes the scaling expo- 
nents to take smaller values corresponding to larger anomalous corrections. 
Again, this situation is reminiscent of the case of the boundary layer. In 
fact, the measured values in our simulation are remarkably similar to those 
reported in the atmospheric boundary layer [|l^ or in a turbulent boundary 
layer They are in between the values two different values measured 

by Toschi (|^) in and above the logarithmic layer in numerical DNS of a 
channel flow. A summary of these results is given in Table When a tur- 
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Figure 8: Same caption as in Fig. |^ for £ = 27r/64 

bulent viscosity is added to the RDT simulation, the intermittent wings are 
less pronounced in the PDF's and the anomalous correction decrease (Table 
|l|), becoming similar to those observed in the DNS. This agrees with the 
picture in which the anomalous corrections are determined by the non-local 
interactions, while the local interactions act to restore the classical Gaussian 
(Kolmogorov-like) behavior. Obviously, the shape of the turbulent viscosity 
also influences the intermittency correction: for the transverse case, where 
there is no asymmetry of the PDF's, both the constant turbulent viscosity 
and the RNG turbulent viscosity provide intermittency corrections which are 
of the same level as the DNS. This is quite remarkable, since they include 
only one adjustable parameter, tuned as to conserve the total energy. For 
the longitudinal case, where an asymmetry is present, the two prescription 
give noticeably different result: as one goes towards lower scales, and as the 
asymmetry becomes larger between the positive and the negative increments, 
the PDF's computed with RDT and constant turbulent viscosity display tails 
which are very close to that of DNS, while the PDF's of the RDT with tur- 
bulent RNG viscosity have a tendency towards a symmetrical shape, thereby 
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Figure 9: Same caption as in Fig. |^for £ = 27r/4 

failing to reproduce the DNS behavior. This difference of behavior between 
the RNG and constant turbulent viscosity will be further investigated in 
Section 2. 

1.4 Comparison of the "local" experiment with DNS 

Given the comparison of the DNS and the RDT ("non-local") simulation, 
one could argue that the increase of the intermittency in RDT is mostly 
due to the increased mean intensity of the small scales (which is seen on the 
energy spectrum plot). A similar increase of the small-scale intensity could be 
produced by other means which have nothing to do with non-locality, e.g. by 
reducing viscosity in DNS. Will there be stronger intermittency in all of such 
cases too? In order to prove that it is not the case, we perform a simulation 
where, as the opposite of the RDT one, the non-local interactions at small 
scales were removed from the NS equation and only the local interactions 
were retained. In order to keep the local interactions which involve scales 
close to the cutoff scales, the velocity and the vorticity fields were split into 
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Figure 10: Comparison of the PDF of the transverse velocity increments 
defined by Eq. (|10D for £ = 27r/256 (10^ statistics at t = 0.48 for the ve- 
locity field from the four simulations: DNS (circles), RDT (crosses), RDT 
+ constant viscosity (diamonds) and RDT + viscosity computed from RNG 
(triangles) 

three parts: the large scales , the medium scales near the cutoff, and the 
small scales. This decomposition is defined in Fourier space as follows, 

u(k) = uz,(k) + u„,(k) + u,,(k), (12) 

Uj(k) = UJis{k) + UJms(M) + i^ssi)^) (13) 

where 



Uis(k) 


= u(k) 


for 


k < kJC 




= 


for 


k > kJC 


Ums(k) 


= u(k) 


for 


kf^l C *\ /c *\ C Jvq 




= 


for 


kc/C < k and k > C k^ 


Uss(k) 


= u(k) 


for 


C *\ k 




= 


for 


k ^ C 



(14) 
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Figure 11: Same caption as in Fig. rU for £ = 27r/64 



Using these definition, the equation for the "local" simulation was the fol- 
lowing: 

dtu(k) + P(u-Vcu)(k) 

- [PK-V^..)(k) + P(u,,-V^,.)(k)]^,>,^j = z/Au(k) (15) 

Where P is the projector operator. In our simulation we choose C = 1.2 and 
the same cutoff scale kc = 24. The result of this simulation are compared to 
the equivalent results from the DNS and the RDT simulation. The energy 
spectra are compared in Fig. [TO. This "local" simulation contain more 



energy at small scales than the DNS even the RDT. The bump of energy 
near the cutoff scale kc is due to the fact that the "local" approximation is 
introduced only for scales smaller than kc. Despite the high level of energy at 
small scales, the solution of this "local" simulation is much less intermittent 
than the equivalent field from DNS and RDT. A comparison of the scaling 
exponents is shown Fig. |l^ and table |I| for both the longitudinal and trans- 



verse velocity increments. These results confirm the idea that intermittency 
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Figure 12: Same caption as in Fig. W for i = 27t/4 



is caused by the non-local interactions and not just by mere presence of the 
small scales. 
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Figure 13: Structure functions computed with the PDF of the longitudinal 
velocity increments from DNS (circles) and the RDT (diamonds) for the first 
(sohd hne), second (dash-dotted line) and third (dash-dot-dot line) moments 
(for convenience of presentation, the RDT structure functions have been 
multiplied by 10) 
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Figure 14: Same plot as in the Fig. O for the transverse velocity increments 



Table 2: scaling exponents of the velocity structure functions from : at- 
mospheric turbulence at 10000 < Rx < 15000 (Dhruva) [Q, channel flow 
(Toschi) near the wall (20 < < 50) and far from the wall {y~^ > 100) 
at Re = 3000 and boundary layer at Rs = 32000 (Zubair) [|5|, |54 





longitudinal 


transverse 


Order 


Dhruva 


Zubair 


Toschi 20<y+ < 50 


Toschi y+ > 100 


Dhruva 


1 


0.366 




0.44 


0.37 


0.359 


2 


0.700 


0.70 


0.77 


0.70 


0.680 


3 


1.000 


1.00 


1.00 


1.00 


0.960 


4 


1.266 


1.20 


1.17 


1.28 


1.200 


5 


1.493 


1.52 


1.31 


1.54 


1.402 


6 


1.692 


1.62 


1.44 


1.78 


1.567 


7 




1.96 


1.55 


2.00 
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Figure 15: Comparison of the scaling exponents computed from the DNS 
(diamonds) and the RDT (circles) statistics. The longitudinal exponents are 
plotted with empty symbols and the transverse exponents with filled symbols. 
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Figure 17: Comparison of the scaling exponents computed from the DNS (di- 
amonds), RDT (circles) and the "local" simulation (squares) statistics. The 
longitudinal exponents are plotted with empty symbols and the transverse 
exponents with filled symbols. 



27 



2 Qualitative explanation of the intermit- 
tency 

Our results can be used to get a qualitative understanding of the intermit- 
tency via the scale behavior of both the PDF of the velocity increments and 
its moments. For this, we are going to build a toy model of turbulence, 
mimicking the small scale non-local dynamics. In spirit, this amounts to an 
"anti-shell" model of turbulence, because, here, we retain only interactions 
between distant wave-numbers, by contrast with the ordinary shell model 
23, which theoretically only retains local interactions. In this con- 



text, it is interesting to note that the ordinary shell model requires a certain 
degree of non-locality between modes so as to generate intermittency: it can 
indeed be proven that the intermittency correction disappear when the sep- 
aration between two consecutive shell tends to zero fl^. Another known 



pitfall of the shell model is its incapacity to describe the observed skewness 
(asymmetry) generation along the scale of the PDF of the longitudinal in- 
crements. This is annoying, since this skewness is directly related to the 
non-zero value of the third order moment, and, hence, to the essence of the 
Kolmogorov cascade picture via the 4/5 law. Finally, the original shell model 
is very crude, since there is no spatial structure (everything is described by 
Fourier modes). We now show how elaborate a cleaner toy model of turbu- 
lence using localized wave-packets, leading to a description of the small-scale 
statistics in term of Langevin processes subject to coupled multiplicative and 
additive noise. 



2.1 The Langevin model of turbulence 

Our numerical simulations showed that both the energy spectrum (and de- 
cay) and the intermittency quantities are well reproduces by a model in which 
only the non-local interactions are left in the small-scale equations whereas 
the local interactions of small scales among themselves are replaced by a 
turbulent viscosity term. Such a model is described by equations (|^) with v 
replaced by the turbulent viscosity coefficient Vt in the small-scale equations, 

dtUi + dj{UiUj) + dj{uiUj) = -diP + ut^ui + di, (16) 

djUj = 0, 
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and (Ji is given by (0). We are interested in the contribution of non-local 
interactions to the statistics of the non-Gaussian small scales. For this, we 
assume that the large scale L quantities ( U and its derivative, and cTj) are 
fixed external processes, with prescribed statistics (to be defined later), and 
derive an equation for the small scale 1 velocity field u', by taking into account 
the scale separation Ij L = e <^\. For this, we decompose the velocity field 
into localized wave-packets via a Gabor transform (GT) (see 



M(x,k,t) =1 (7(e,|x-x'|)e*('^-^')u(x',t)dx', (17) 

where (7 is a function which decreases rapidly at infinity and 1 ^ <^ e. 
Note that the GT of u is a natural quantity for the description of the velocity 
increments because of the following relation 

u (x + 1) - u (x - 1) = ^ y" e-^-^lm u(x, k) dk. (18) 

Thus, velocity increments are related to GT via the Fourier transform, and 
all information about the /-dependence is contained in the GT dependence 
on k (the main contribution to the above integral comes from k ~ 2t^/1.). On 
purely dimensional ground, we see that u ~ k~'^Su, where d is the dimension. 
Therefore, in the sequel, we shall identify k'^Imu with the velocity increment. 
Applying GT to (H) we have (see ||38| for details): 

Dtu = u ■ ^ + (Ti. - utk^u, (19) 

where ^ and a± are random processes, given by 

i = v(2^U.k-u), 

ax = a--^(k-a), (20) 



and 



k 



Dt = 9t + i-V + k-Vk, 

i = V = VkH, (21) 

k = -V(k ■ U) = -VH, (22) 

H = U-k. (23) 
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Because the large-scale dynamics is local in fc-space, it is only weakly af- 
fected by the small scales and the quantities ^ and cr in equation ([T9|) can 
be considered as a given noise. Also, because the equation is linear in it, we 
immediately see that k'^u will also satisfy an equation similar to ([19|) subject 
to a straightforward modification of the force definitions. Before elaborating 
more on (|19]), it is convenient to study in closer details the physical param- 
eters of this equation. 



2.2 The noises 

In equation (0), the noise ^ and a± appearing as the projector of two quan- 
tities, one related to the velocity derivative tensor, and another related to 
the Gabor transform of the energy transfer from large to small scales. In the 
sequel, we present a statistical study of these two noises in the physical space 
(i.e. after inverse Gabor transforming a± ). We will also consider the Fourier 
spectra of the corresponding two-point correlations. Let us choose k to be 
along one of the coordinate axis (without loss of generality because of the 
isotropy). Then the components of C, coincide (up to the sign) with the cor- 
responding components of the velocity derivative tensor, and we will use this 
fact in the rest of this section. The velocity derivative tensor has been studied 
in the literature e.g. in 0] in terms of correlations between the directions 



of the anti-symmetrical part of the tensor (the vorticity) and the symmetrical 
part (the strain). Over long time, the vorticity appears to be aligned with 
the direction of largest stretching. Other studies focused on the PDF of the 



modulus of one component. For example, Marcq and Naert observe that 
the derivative has a highly non-Gaussian distribution, but with a correlation 
function which decays rapidly, and can be approximated by a delta function 
at scales large compared to the dissipative scale. In the present case, we ob- 
serve different features. Because of the isotropy, we can concentrate only on 
two quantities, say and ,^12. Fig. shows the equal-position, time corre- 
lation Cii(t — to) =< ^ii(t)^ii(to) > and ajj(t — to) =< o'i(t)ai(to) > as a func- 
tion of t — to- Note that these quantities are normalized to 1 at t = to in this 
figure. Firstly, we see that Cn approximately coincide with C12 and au co- 
incides with which is a good indicator of isotropy (without normalization 
there would be C12 = —3 Cn). Secondly, we see that the correlations Cu and 
C12 decay to zero over a time scale which is of the order of few turnover times 
r (r = 0.19). On the other hand, the correlations of a decay much faster. 
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<CT,(t),a,(t„)>: 

<CT,(t),a,(t„)>: 



0.5 1.0 1.5 2.0 2.5 



Figure 18: The normalized auto-correlation in time for the two force compo- 
nents a and ^ (the turnover time r is equal to 0.19) 

over a time of the order of r/2. Fig. |19| displays the Fourier transforms of the 
equal-time two-point correlations Dxiix—x^) =< C,ii{x, y, z, t)^ii{xo, y, z, t) > 
and aii{x — Xo) =< crj(x, t)(Tj(xo, t) >. One sees that all correlations are very 
weak beyond 2kc- The correlation D12 appears to be the largest at large 
scales, but it decays more rapidly than the other correlations. We have also 
investigated the cross-correlations between the noises. The equal position 
cross-correlation are displayed in Fig. ^ The correlation is rather weak, 
but there is a tendency for cTj to be correlated with over a time scale of the 
order of r, while it is anti-correlated with the other component of the tensor, 



over such a time scale. The equal time correlations are shown in Fig. ^ 
Notice that the cross correlation are one order of magnitude weaker than the 
direct correlations. The cross-correlation involving ^12 are essentially zero, 
while the correlations involving display a first overall decay up to A; = kc, 
followed by an extra bump up to the end of the inertial range {k = 40). In 
the next Section, we will show that this feature is actually related to the 
energy cascade. 
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<e,z(x),e„(x„)> 




Figure 19: Fourier transform of the space auto-correlation for the two com- 
ponents of the force a and ^. Coordinates y and z are fixed. 

One may also note that the two noises are spatially very intermittent. In 



Fig. we show iso-surfaces of the modulus of the noises, corresponding 
to 3.5 times the RMS value. For comparison, the same plot is made for 
the vorticity. One observes well defined patches of cr which are strongly 
correlated with areas of strong vorticity. In the case of ^, the patches are 
much more space filling. The longitudinal component is characterized by 
smaller-scale structures than the transverse component ^12. 

To obtain an indication about the scale variation of the statistical proper- 
ties of the noises, we also computed the PDF's of the noise spatial increments 
ban = a{x + i) — cr(x) and S^ij^, = ^ij{x + i) — iij{x). Note that the first 
of these quantities is directly related with the Gabor transform of the addi- 
tive noise, whereas the second one contains some useful information about 
the time correlations via the Taylor hypothesis (which is valid locally be- 
cause the large-scale velocity is typically greater than the small-scale one). 



Fig. |23| shows the results of the longitudinal and transverse increments for 



the first component of the additive noise Oi and the equivalent results for 
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Figure 20: The cross-correlation in time between the forces a and ^ (t = 0.19) 



the component of the muhiphcative multiphcative noise are shown on 



Fig. |2J. The PDF are displayed for i = 27r/256 and 27r/4. One observes 



wide, quasi-algebraic tails for the additive noise, similar to those observed 
for the PDF's of the velocity derivatives. The PDF of are much closer to 
Gaussian statistic. 



2.3 The turbulent viscosity 

In the previous Section, we have discussed the influence of two prescription for 
the turbulent viscosity, one based on the RNG, one taken simply as constant. 
In the sequel, we shall use the simple formula: 

where and B are a constant and u = k'^ Imu is the velocity increment 
over a distance 1/k (hereafter we drop 5 in 5u). When 5 = 0, this for- 
mula provides the constant turbulent viscosity. When i/q = 0, it provides a 
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Figure 21: Fourier transform of the space cross correlation between the force 
a and the force ^. Coordinates y and z are fixed. 

dimensional analog of the RNG viscosity, and tends to zero as k tends to 
infinity. 



2.4 Statistical properties of the velocity increments 

We are now going to derive qualitative results by adopting two complemen- 
tary points of view: in a first one, we will study the statistical properties of 
the velocity increments in the frame of reference moving together with the 
wave packets in (fc, x) space. This corresponds to a Lagrangian description in 
the scale space. In the second approach, we replace time with its expression 
in terms of k, as it would follow from the ray equation (p2D. This will give 
as an equation at a fixed k which corresponds to an Eulerian description. As 
a further simplification, we shall leave for further study the possible corre- 
lation between longitudinal and transverse velocity increments described for 
example in |^ and consider a one-dimensional version of (|1^), treating the 
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Figure 22: Isovalue (3.5 times the RMS) at to = of a) the absolute value 
of vorticity, b) the corresponding additive force (|(t|) c) .^n = dUx/dx and 
d) 62 = dUjdy 

quantity u = k Im XL cLS cL "velocity increment" over the distance / = 27r/A;, 

Dtu = + a± — utk'^u, (25) 
k = -k^ (26) 

Here, we assumed that the forcing to be symmetric such that it does not 
produce any Reu. This toy model can also be viewed as a passive scalar in 
a compressible ID flow. Artificial introduction of compressibility is aimed at 
modeling the RDT stretching effect which appears only in the higher number 
of dimensions for incompressible fluids. 

Study of the noises in Section 2.2 revealed their rich and complex behav- 
ior. As a first simplifying step, we disregard these complexities and use the 
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Figure 23: PDF of the increments + z) — (Ti(x, 2;) for ix = 27r/256 
(circles) and ii = 27r/4 (crosses) and cri^x.y + £21-2) — o-i{x,y,z) for 
£2 = 27r/256 (triangles) and £2 = 27r/4 (diamonds) (the dotted line corre- 
sponds to Gaussian statistics). 
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Figure 24: PDF of the increments ^ii{x+£i, y, z)—^ii{x, y, z) for £i = 27r/256 
(circles) and £i = 27r/4 (crosses) and S,ii{x,y + ^2,2;) — ^ii{x, y,z) for 
£2 = 27r/256 (triangles) and £2 = 27r/4 (diamonds) (the dotted line corre- 
sponds to Gaussian statistics). 

Gaussian, delta correlated approximation, as will be done in the next two 
sub-sections. Given a rather short time correlation of a, our delta approxi- 
mation is rather safe. The delta approximation for ^ is more debatable, and 
the performance of such a model should be further examined in future. Also, 
the Gaussian hypothesis is obviously only valid at large scale, and for .^n. 
Therefore, the generalization of ours results for non-Gaussian noises would 
be very interesting, and is the subject of an ongoing research. In the sequel, 
we consider the function a, D and A as free parameters. 

2.4.1 The Lagrangian description 

In the frame of reference moving with the wave packets in (A;, x) space, the 
l.h.s. of ( p5D becomes simply the time derivative. On the other hand, k has 
to be replaced in terms of its initial value /cq and time everywhere including 
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the noises a and ^. Such a transformation from the laboratory to the moving 
frame can obviously change the statistics of a and ^. In the Lagrangian 
description we will assume that we deal with noises which are Gaussian in 
the moving frame with correlations functions 

< a^{k{t),t)a^{k'{t),t') > = 2a5{t-t'), 
<mt),mk'{t),t')> = 2D6it-t'), 
<i{k{t),t)a^{k\t),t')> = 2\S{t-t'), (27) 

where coefficients a, D and A depend on the scale via ko. With these noises, 
(p5D becomes a Langevin equation for the velocity increments, where ^ is 
a multiplicative noise, a± is an additive noise, and u^k'^u the ( non-linear) 
friction. The multiplicative noise is produced by interaction of two small 
scales with one large scale whereas the additive noise is due to a merger 
of two large scales with into one small scale (therefore, the later acts at 
the largest among the small scales). For Gaussian, delta correlated noises, 
this Langevin equation leads to a Fokker-Planck equation for the probability 
distribution Pk{u,t) of the velocity increment u, 

dtPk = du (^Uik'^uPk) 

+ Dkdu {uduuPk) - Xkdu (uduPk) 

- XkdliuPk) + akdlPk, (28) 

where we have taken into account the fact that, due to homogeneity, and a 
have a zero mean. Here, we dropped the the subscript in /cq and dependence 
of the all involved quantities on the scale is simply marked by the subscript k 
(the scale dependence a, D and A is still unspecified). The stationary solution 
of dD is: 

-vtk'^y - Dy + X 
/o Dy^ -2Xy + a 

where Ck is a normalization constant. The integral appearing in ( p9|) can be 
explicitly computed in two regimes: in the first one, for -u << kuo, we have 
(z/( = z/q), and we simply get 

C 

^'^""^ = (Z^«2-2A« + a)^/^+-^^/(^^)' 



P,(n) = C,exp/ dy, (29) 
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The range of u for which the PDF follows this algebraic law decreases with 
increasing scales. It is the largest (and hence it is best observed) at the 
dissipative scale, where the velocity increments are equivalent to velocity 
derivative or to vorticity. Several remarks are in order about this expression. 
First, notice that the distribution is regularized around m = by the presence 
of the parameter a, but then displays algebraic tails. These are well known 
features of random multiplicative process with additive noise (see e.g p5| ). 
The occurrence of algebraic tails in vorticity PDF's has been noted before by 



2J, ^ in the context of 2D turbulence. However, processes with algebraic 
tails are characterized by divergent moments. These divergences can be re- 
moved by taking into account finite size effects, like physical upper bounds 
on the value of the process (see, e.g. [^] for discussion and references) which 
introduce a cut-off in the probability distribution. This effect is automat- 
ically taken into account in our simple model, via the turbulent viscosity, 
which prevents unbounded growth of velocity fluctuations and introduces an 
exponential cut-off. 

Indeed, in the regime u» kuQ, we see that: 



dlnP —u\u\ 



(31) 



du Du^ — 2Xu + a 
This mean that at large u the PDF decays like an exponential, or even faster 



if D = (see below). This feature had been noted by Min et al and 
finds here a detailed explanation. 

Another important observation is that the PDF's have an intrinsic skew- 
ness, which can be traced back to the non-zero value of A, i.e. to the correla- 
tion between the multiplicative and the additive noise. The physical picture 
associated with this correlation is related to the correlation between vortic- 
ity (present in the large scale strain tensor) and stretching (associated with 
the term U"VU, present in a), which is the motor of the energy cascade [|I|]. 
The importance of the additive noise in the skewness generation has been 



stressed before ||T6[. We find here its detailed explanation. Note also that 
the trends towards Gaussian large scale behavior of the velocity increments 
can be easily accounted for if the multiplicative noise tends to zero at large 
scale {D, A — > 0). In such case, the process becomes purely additive, and the 
limiting PDF is Gaussian for u independent turbulent viscosity, or can also 
fall faster than a Gaussian (like exp(— M'^),for Ut (x u). Such a supra Gaussian 
behavior has been noted before In between the dissipative and the largest 



39 



scale, the transition operates via PDF's looking like stretched exponential. 

This qualitative feature can be tested by comparison with the numerical 
PDF's. Our model, predicts that dlnP/du should behave like the ratio: 



d\nP -u\l iniiuY + ml- rUiU + ma 



' ^ " (32) 



du rri/^u^ — 2m^u + 

Without loss of generality, we can f actor ize out the parameter m^. The fit 
Therefore only contains four free parameters, which can be easily related to 
the physical parameters of the problem. We have computed this derivative 
for the PDF of longitudinal and transverse velocity increments at various 
scale, and performed the four parameters fit. Examples are shown on Fig. 



25| . Observe the good quality of the fit, but we stress that there is a rather 
large uncertainty in the determination of the parameter of the fits, which 
sometimes cannot be determined better than up to a factor of two by our 
fit procedure (a standard least square fit). The scale dependence of the 
coefficients of the fit, is shown on Fig. p6| and ETl Note that since transverse 



velocity increments are symmetrical by construction, we have set ma = 
in the fit. Note that in the longitudinal case, ma is smaller than the other 
coefficients by two orders of magnitude. This features the smallness of the 
skewness, and the weakness of the correlation between the multiplicative 
and the longitudinal noise (see Section 2.1). The parameters have a rough 
power law behavior (see Fig. |26|and^^. Theoretically, one expects the ratio 



mi/ 1712 to behave like l/(z/o/c), if ( p^ ) holds. The power- law fits of Fig. |26 
and ^ provides a dependence of k^'^^ for the longitudinal case, and 
for the transverse case, corresponding to a scale dependence of uq given by 
k^-^^ and 

2.4.2 The Eulerian description 

We now consider (^), again in the moving with the wave packet frame, but 
now we change the independent variable from t to k = k{t) which satisfies 
(HD. We get 

-Vu H — , (33) 



d\iik i i 

where "P is a number, accounting for the projection operator (a way to mimic 
higher- dimensional effects in our ID toy model). Eq. (|33| ) is again a Langevin 



40 



n 



Figure 25: Fit of dlnP/du for the transverse increments, at 1/k = 2 
(squares), 1/k = 32 (diamonds), 1/k = 42 (circles) and 1/k = 52 (trian- 
gles). The fits via formula ^ are given by lines. Note the good quality of 
the fit. Similar fits for other scale separations, and in the longitudinal case 
were obtained. 

equation for velocity increments in the scale space, with multiplicative and 
additive noises which is now expresses in an Eulerian form. Similar Langevin 
processes have been proposed before to explain the scale dependence of ve- 
locity increments [^, but without additive noise |0 . 

The noises in this Langevin equation are different from the noises ap- 
pearing in the Lagrangian representation and they would have a complicated 
statistics if we assumed that a and ^ were Gaussian in the Lagrangian repre- 
sentation. However, we can simply assume here that the noises are Gaussian 
and delta-correlated in the Eulerian representation (which is different from 
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the assumption of the previous section) and re-define a, D and A as 

< t{k))C^{k, t{k))a^{k', t{k'))C\k', t{k')) > = 2a 6{k - k'), 

<r\k,t{k))r\k',t{k'))> = 2D6{k-k'), 
<r\k,t{k))a^{k',t{k'))r\k',t{k'))> = 2X6{k-k'),{M) 

This allows one to derive the Fokker-Planck equation corresponding to (|33|) , 
kdkP{u,k) = duiVuP) 

+ Ddu {h'tk'^udu vtk'^uP ) - Xdu {I'tk'^uduP^ 

- Xdl (vtk'^uP) + adlP. (35) 



We may use (|35D to derive an equation for the moments, by multiplication 
by and integration over u. With the shape of the turbulent viscosity given 
by (H, we get: 



kdk < > = -C(n) < > 



+ nDB^e < > -\n{2n - l^k^ < > 



- 2n^\B^ < > +an{n - 1) < u""^ >, (36) 

where ({n) = nV — n^Dk'^v'^ is the zero-mode scaling exponent. For n = 1 
and taking into account the constraints that < u >= (homogeneity), one 
get a sort of generalized Karman-Horwath equation: 

As in the Lagrangian case, this means that skewness (related to non-zero < 
> is generated through non-zero value of A, i.e. through correlations of the 
multiplicative and the additive noise. However, due to the turbulent viscosity, 
we cannot explicitly solve the hierarchy of equation. In many homogeneous 
turbulent fiows, however, the skewness (proportional to A) is quite small, and 
moment of order 2n+l are generally negligible in front of moments of order 
2n + 2. For even moments, this remark suggest that to first order in the 
skewness, and for 2n > 1 the dynamics is simply given by: 

kdk < M^" > = -C(2n) < M^" > +2nDB^e < u^''+^ > 

+ a2n{2n - 1) < w^^-a > ^o{\^). (38) 
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Note that D/a given by the parameter m4 in our fit, and is such that 
Dk"^ /a increases with k (see Figs. ^ and |2^). Therefore, at small scales. 




Figure 26: Coefficient of the fit of dhi{P)/du with (|32D as a function oil/k 
for longitudinal velocity increments: mi (squares), m2 (empty diamonds), 
(circles) and m4 (filled diamonds). The lines are the power-law fits: 
(long-dashed line); (short-dashed line); (dotted line) and 

(dash-dot line). 

the dominant balance is: 

kdk < M^" >= 2nDB^k^ < > . (39) 

The solution is < m^" >oc fc"^". This is the usual "regular" scaling in the 
dissipative zone. For larger scales, DB^k^ -C a, and if a varies like a power 
law, the general solution of (|38| ) is a sum of power-laws: 

n 
p=0 
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Figure 27: Coefficient of the fit of d\'n.{P)/du with (|32D as a function of 1/k 
for transverse velocity increments: mi (squares), m2 (diamonds), and m4 
(circles). The lines are the power-law fits: k~^'^^ (long-dashed line); 
(line); and (dotted line). 



This solution illustrates the famous mechanism of "zero- mode intermittency" 



[ P^ . Here, the zero mode is the solution of the homogeneous part of eq. (pq), 
i.e. a power-law of exponent —({2n). Without the zeroth mode, responsible 
for the ffist n — 1 scaling laws, the 2n-th moment will scale in general like 
a (A;)" (provided one assume that this dominates the other term), i.e. will be 
related to the turbulent forcing. This is the standard Kolmogorov picture. 
When the zeroth mode is taken into account, the moment now includes new 
power laws, whose exponent is independent of the external forcing, and which 
can be dominant in the inertial range, thereby causing anomalous scaling. In 
the present case, the scaling exponent is quadratic in n, and reflects the log- 
normal statistics induced by the Gaussian multiplicative noise, in agreement 
with the latest wavelet analysis of Arneodo et al M. Note also that the 
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competition between the zeroth mode scaling and the scahng due to external 
forcing forbids the moments to scale like a power-law, thereby generating a 
breaking of the scale symmetry. 

For odd moments, we cannot perform any rigorous expansion because all 
the terms of the equation are of order A. For low order moments, however, 
the computation of < u^^^/vt > mainly involve velocity increments close to 
the center of the distribution, for which z/j ^ uq. So, for low order, it is 
tempting to approximate the equation for the odd-order moments by 

kdk < ^2"+^ >^ -C{2n + 1) < > -2n^\B^ < u^''+^ > u-\ (41) 

This approximation is only valid in the inertial range, where the last term of 
( |36D can be neglected. An immediate consequence of this loose approximation 
is that C{2n + 1) = C{2n + 2) in the inertial range. Odd moments (without 
absolute values) are very difficult to measure because of cancelation effects 
which introduce a lot of noise. In our case, due to our limited inertial range, 
we cannot compute these exponent with a sufficient degree of accuracy. A 
careful investigation performed in a high Reynolds number boundary layer 
however seem to be in agreement with our prediction, as is shown in Fig. 
|28| : it is striking to observe that ({5) ~ C(6), C(7) ~ C{^)y stc making the 
curve look as if odd and even scaling exponents are organized on a separate 
curve [Q. A second independent experimental check of our prediction (p]) 
is that (C(2n + 2) - C(2n + 1)) < u^''+^ > / < m^^+s > should scale, forn > 3 
like n^. Fig. ^ shows that this is indeed the case. 

3 Discussion 

In this paper, we have shown that non-local interaction are responsible for 
intermittency corrections in the statistical behavior of 3D turbulence. Re- 
moval of the local interaction in numerical simulations leads to a substantial 
increase in the number of the tornado-like intense vortex filaments and to 
a stronger anomalous corrections in the higher cumulants of the velocity in- 
crements. It is also accompanied by a modification of the energy transfer 
in the inertial range, tending to create a flatter energy spectrum. The in- 
termittency corrections and the spectra are close to that observed in high 
Reynolds number boundary layer, suggesting that the non-local interactions 
prevail in this geometry. This could be explained by presence of the mean 
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Figure 28: Exponents of the structure function in a high Reynolds number 
boundary layer |T^. Note the tendency for ({2n + 1) = ({2n + 2) for n > 3. 



flow, which geometrically favors non-local triads in the Fourier transform of 
the non-linear interactions. We showed that replacing the removed local in- 
teractions by a simple turbulent viscosity term allows to restore the correct 
intermittency and the energy characteristics. Our results agree with the be- 
lief that intermittency is related to thin vortices amplified by the external 
large scale strain similarly to the classical Burgers vortex solution. Local in- 
teractions can be viewed at mutual interaction of these thin intense vortices 
which result in their destruction, which is also in agreement with our results. 

To prove that the enhanced intermittency is not simply the result of 
the stronger small-scale observed in the RDT simulation, we performed yet 
another numerical experiment in which the non-local interactions were ne- 
glected and only the local ones retained. This resulted in even higher (than 
in RDT) level of the small scales but it exhibited much less intermittency, 
which confirms our view that non-locality is crucial for generation of the 
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Figure 29: (C(2n + 2) — C{2n + l))A2n+i/A2n+2 as a function of n in a 
high Reynolds number boundary layer. Here An is the prefactor of the 
(non-dimensional) structure function of order n. The line is n^, the pre- 
diction of our toy model. 

intermittent structures. 

The result that the net effect of the local interactions is to destroy the 
intermittent structures is at odds with a very common belief that the inter- 
mittency is due to the vortex reconnection process which takes a form of a 
finite time vorticity blow-up. Indeed, the later is a strongly nonlinear pro- 
cess in which the local vortex- vortex interactions are important. However, 
this process seems to be dominated by another local processes the net re- 
sult of which is to destroy the high- vorticity structures rather than to create 
them. It would be premature to claim, however, that the same is true at any 
arbitrarily high Reynolds number. 
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Our numerical approach sets severe limitations to the value of the Reynolds 
number we are able to explore. In this context, it is interesting to point out 
that preliminary tests regarding the importance of non-local interactions have 
been conducted on a velocity field coming from a very large Reynolds num- 
ber boundary layer |]lOl- Even though the test is not complete (the probes 
only permit the accurate measurement of special components of the veloc- 
ity field), it tend to suggest that non-local interactions dominate the local 
interactions by several orders of magnitude. Our results would also explain 
the findings of the Lyon team who found that when probing fluid area 
closer and closer to a large external vortices, or to a wall boundary, one could 
measure energy spectra moving from a k~^/^ law towards a spectra, while 
anomalous corrections in scaling exponents would become more pronounced. 
At the light of our study, this could be simply explain by a trend towards 
more non-local dynamics via the mean-shear effects at the wall. 

Based on the conclusions of our numerical study, we developed a new 
toy model of turbulence to study the intermittency. It has the form of a 
Langevin equation for the velocity increments with coupled multiplicative 
and additive noise. We showed how this toy model could be used to under- 
stand qualitatively certain observed features of intermittency and anomalous 
scaling laws. Among other things, we showed how the coupling between the 
two forces is related to the skewness of the distribution, and how algebraic 
and stretched exponential naturally arise from the competition between the 
multiphcative and the additive noise. We tested our qualitative predictions 
with experimental and numerical data, and found good general agreement. 
To be able to turn our toy model into a tool for "quantitative" study of the 
intermittency, several developments are needed. The first one is to consider 
the multi-dimensional version of our toy model, to be able to couple longi- 
tudinal and transverse increments. The scale dependence of the turbulent 
viscosity and of the forcing needs to be further investigated, possibly using 
tools borrowed from the Renormalization Group theory (see e.g. ||^). Also, 
the non-Gaussianity of the noises could be taken into account. 

In 1994, Kraichnan proposed an analytically solvable, new toy model 
for the passive scalar, which provided a substantial increase of our under- 
standing of the passive scalar intermittency. Our model, built using the 
non-local hypothesis, is a direct heir of this philosophy in that, as in passive 
scalars, all important intermittency effects are produced via a linear dynam- 
ics. The nonlinear (local) scale interactions are important too because they 
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are the main carrier of the energy cascade, but it is only their mean effect 
and not statistical details that are essential, 
the present paper. 
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